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Abstract 

We study the three-dimensional U(l)+Higgs theory (Ginzburg-Landau model) as an 
effective theory for finite temperature phase transitions from the 1 K scale of super- 
conductivity to the relativistic scales of scalar electrodynamics. The relations between 
the parameters of the physical theory and the parameters of the 3d effective theory 
are given. The 3d theory as such is studied with lattice Monte Carlo techniques. 
The phase diagram, the characteristics of the transition in the first order regime, and 
scalar and vector correlation lengths are determined. We find that even rather deep 
in the first order regime, the transition is weaker than indicated by 2-loop perturba- 
tion theory. Topological effects caused by the compact formulation are studied, and 
it is demonstrated that they vanish in the continuum limit. In particular, the photon 
mass (inverse correlation length) is observed to be zero within statistical errors in the 
symmetric phase, thus constituting an effective order parameter. 
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1 Introduction 



Finite temperature phenomena are very important for cosmology and for the physics 
of uhrarelativistic heavy ion colhsions. It is thus indispensable to have efficient and 
accurate methods for computing, e.g., the partition function for a given theory. This 
paper is devoted to solving this problem for the U(l) gauge+Higgs theory, or scalar 
electrodynamics, or the Ginzburg-Landau model. This theory in three dimensions is 
the effective theory of superconductivity [|1| and of liquid crystals in certain regimes, 
and it is also an interesting theoretical laboratory for studying the thermodynamics of 
and topological defect formation in relativistic field theories. 

A finite temperature system singles out a specific Lorentz frame, the rest frame of 
the matter. In non-relativistic condensed matter physics this automatically leads to a 
three-dimensional (3d) — or in some cases even lower dimensional — formulation of 
thermodynamical computations. The Hamiltonian is given and the partition function 
is computed in 3d. For realistic systems the full computation is usually too complicated 
and one has to replace the original theory by an effective one focussing on the essential 
degrees of freedom. For instance, in superconductivity the full quantum theory of 
electrons in an ionic lattice is first replaced by the BCS theory, which then for a class 
of phenomena can be replaced by the very simple Ginzburg-Landau model 

The situation is different in finite temperature relativistic field theories. The general 
first-principles Lorentz and gauge invariant formulation is necessarily four-dimensional 
(4d). However, even here, for a class of theories and phenomena, the full 4d theory can 
be replaced by effective 3d theories of various degrees of simplicity. The original idea [3- 
5] is rather old and has recently been applied to a large class of relativistic field theories, 
such as QCD [6-9], the SU(2)+fundamental Higgs theory or the Standard Model [10- 
13], the Minimal Supersymmetric Standard Model [14-17], and the SU(5)+adjoint 24- 
plet Higgs theory |18|. The corresponding 3d effective theories have been numerically 
studied for QCD in [19-21,9], for SU(2)+fundamental Higgs in [10,22-27], and for 
SU(2)xU(l) + fundamental Higgs in The present case of the U(l)+Higgs theory, 
for which dimensional reduction has been considered in [12,29-31] has been numerically 
studied in [32-38]. 

The purpose of this paper is to study the U(l)+Higgs theory in the effective theory 
approach, emphasizing the two essential but totally separate aspects of the problem: 

• The numerical study of the 3d effective theory as such: the determination of the 
phase diagram, of the characteristics of the phase transitions found, and of the 
correlation lengths of various gauge invariant operators. We have presented our 
first results for the correlation length measurements and for the structure of the 
phase diagram already in and here we concentrate on the characteristics of 



the phase transition in the 1st order regime and on explaining the details of the 
correlation length measurements. Earlier numerical work has appeared in [32-37]. 
The full solution of the problem requires this numerical work, although analytic 
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methods [39-44] give important insights, as well. It is also quite interesting to 
compare the U(l)+Higgs theory with several related theories with a massless 
photon in the tree-level spectrum: SU(2)xU(l) + fundamental Higgs theory 
[0, SU(2)+adjoint Higgs theory ^ § and the pure U(l) theory [45-47]. The 
topological mass created for the photons in the compact formulation [48-51] will 
play an important role in this context. 

The analytic derivation of the relation between the full physical theory and the 
effective theory. Here many different physical theories can be mapped onto one 
effective theory. The U(l)+Higgs theory is particularly interesting in that the 
full theory can be, e.g., superconductivity at T ~ 1 Kelvin or hot scalar+fermion 
plasma at relativistic temperatures. Of course, the latter case can also be ap- 
proached directly from the 4d viewpoint [52, 53[ . 



The results for the first aspect — the numerical study of the 3d theory as such — 
are contained in Sections 2-6 of this paper. We formulate the problem in continuum 
in Sec. 2, make some perturbative computations in Sec. 3, discretize the theory in 
Sec. 4, discuss the photon in the discretized theory and in some other 3d theories in 
Sec. 5, and present our main lattice results in Sec. 6. Since the second aspect of the 
problem, the derivation of the 3d theory, completely factorizes from the first aspect, 
we discuss the derivation separately in Appendix A for the case of superconductivity 
and in Appendix B for the hot scalar -l-fermion plasma. We conclude in Sec. 7. 



2 3d U(l)+Higgs theory in continuum 

The 3d U(l)+Higgs theory is a locally gauge invariant continuum field theory defined 
by the functional integral 



Z = j VAiVcf) exp[-5(Ai 
■1 



where Fi,, 



S 



diAj djAi 



F"^- + ID- 



+ ml(f)*(j) + A3 



(2.1) 
(2.2) 



and Di = di + ie^Ai. The parameters ms. 



e|, A3 of the 



Lagrangian have the dimension GeV and the fields have dimension GeV^^^. Since 
the theory in eq. (2^) is a continuum field theory, one has to carry out ultraviolet 
renormalization. In 3d the couplings e\ and A3 are not renormalised in the ultraviolet, 
but there is a linear 1-loop and a logarithmic 2-loop divergence for the mass parameter 
m|. In the MS dimensional regularization scheme in 3— 2e dimensions, the renormalized 



mass parameter becomes fl^ 



m^(/i) 



Ael + 8A3ei 
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Figure 1: The phase diagram of the G-L theory. The photon mass acts as an effective 
order parameter. 



where Km is a scale independent physical mass parameter of the theory. Instead of it, 
it is more convenient to use 7713(63). Choosing e\ to set the scale, the physics of the 
theory will depend on the two dimensionless ratios 

mliel) A3 
y = ^^. x = -|. (2.4) 

63 63 

We shall later in Appendices A and B relate the parameters of superconductivity and 
finite T scalar electrodynamics to x and y. 

It is now a mathematical and computational problem to determine the physics as- 
sociated with the theory in eq. ( ^.1[ ). As usual all physics lies in expectation values 
of various operators. Since this is a gauge theory, only gauge invariant operators have 
non-vanishing expectation values. The most relevant of these are the local (depending 
only on one point x) operators of lowest dimensionalities (half-odd dimensions could 
be made integer by multiplying by 63): 

• Dim = 1: the J^^ = 0++ scalar 5(x) = 0*(x)0(x). 

• Dim = the vector Vi(x) = Bi = ^eijkFjki^x). 

• Dim = 2: the 1 vector l^i(x) = Im</)*(x)Dj0(x); the l""*" vector S'j(x) = 
Re0*(x)A0(x) = (9i0*(x)0(x); and the 0++ scalar (0*0)^. 

• Dim = 2^: the vector (j)*Bi(f). 

• Dim = 3: The 0^^ scalars FijFij and cffDiDicj); and the 2^+ tensor (f)*[{Di, Dj} — 
21 d 5ijDkDk](t)] etc. 
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• Dim = 3^: The scalar Bidi(f)*(f); and the O""*" scalar Bilm(f)*Di(j); etc. 

• Dim = 4: The 0"*"~ scalar di(j)*(f)lm(f)*Di(j); etc. 

The quantum numbers here refer to 0(3). From these one can further construct, e.g., 
bilocal operators such as correlators of two of the above operators acting at different 
points. 

The first topic of importance is the phase structure of the theory. The critical curve 
y = yc{x) (see Fig. 0) divides the plane into two disjoint regions, the symmetric phase 
ai y > ydx) and the broken phase at y < ydx). The presence of a critical curve is 
signalled by singularities in the free energy 

Z = exp[-Velf{y,x)], (2.5) 

where / is dimensionless. The critical curve is localised by measuring different local or 
bilocal expectation values of the gauge invariant local operators listed above. There is 
no local gauge invariant order parameter, the expectation value of which would vanish 
in one of the phases. Instead, we shall use as an effective order parameter the photon 
mass, which is measured from a bilocal operator. 

For small x, the system has a first-order transition (Fig. p. It then has two bulk 
states of the same free energy, and the broken and symmetric phases coexist at yc- The 
system has one stable and one metastable branch for y-{x) < y < y+{x). Expectation 
values of various gauge invariant scalar operators jump when crossing ydx). The 
following quantities are of particular interest: 

• The jump Afa of the order parameter like quantity £3 = {(f)'^ (l){el)) / between the 
broken and symmetric phases at yc- In perturbation theory, A£s ~ i0^(?/c)/e|, 
where 06 is the location of the broken minimum in, say, the Landau gauge. Note 
that {(j)^ is scale dependent [Q, as given below in eq. ([4.6|). 

• The interface tension = a/{Tef), defined in perturbation theory by 

as = J^^'^^' d{<P/e;)pV{ct>/e,)/el (2.6) 

where V is the perturbatively computed 3d effective potential. On the lattice 
is determined from the fundamental relation F = —pV + a A, which holds in the 
presence of a planar interface separating two coexisting phases. 

In a second order transition, on the other hand, the quantities of interest are the 
different critical exponents. 

The second essential property of the U(l)+Higgs theory is its excitation spectrum. 
In perturbation theory, one can easily see what the free field degrees of freedom are both 
in the symmetric and broken phases (see also Table 1). However, the symmetric phase 
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is in some respects non-perturbative (since the Coulomb potential is logarithmically 
confining in 2+1 dimensions) and thus these states need not be the physical ones. The 
physical states are described by the correlators of the gauge invariant operators listed 
above. For example, the physical mass, or inverse correlation length, of the lowest 
scalar 0"'""'" excitation is obtained from the large distance behaviour of the correlator 
(5'(x)5'(y)). To study vector 1 excitations we shall use the operators Vi(x), V^(x). In 
condensed matter context these masses are usually called "renormalised" masses; we 
prefer to reserve the word renormalisation for the elimination of ultraviolet divergences. 

A particularly interesting state is the photon. It is evidently one of the fundamental 
fields of the action in the symmetric phase; is it in the physical spectrum of the theory? 
In [3S] we have provided non-perturbative evidence for the fact that it indeed is. It 
can thus serve as an effective order parameter and the phase diagram can contain two 
disconnected regions as in Fig. 1. This is a nontrivial fact, since massless states can 
disappear from the physical spectrum or in the lattice formulation of the theory. This 
will be discussed in more detail in Section 5. 



3 Perturbative results for the U(l)+Higgs theory 

Let us consider the region of small x. Then perturbation theory is expected to be 
accurate and it predicts a "fluctuation induced" or "Coleman- Weinberg-type" flrst 
order phase transition. Renormalised perturbation theory induces a dependence on a 
scale /i and in 3d one has to go to 2 loops in order to optimize this scale. The 2-loop 
potential has been discussed in |T^ ^ ^ and optimisation in |T^, and we 



refer there for details (see also eq. (|B.33|) below for the 2-loop potential in the full 4d 
theory with fermions). 

To show the main qualitative features of the perturbative results, consider the 1-loop 
potential in the Landau gauge (see eq. (|B.31|) ), 



^l-loop/ 63 



^?/(/i)02 + ix0^- 



127r 



2^^ + {y + 3x^y/' + {y + x^Y'~2y'^']\. (3.1) 



where cj) = (^/e^ and 



ix 



ix 



In^. 



(3.2) 



The /i dependence of ?/(/i) is of 2-loop order so that it is a higher order effect in eq. ( p.ll) . 
The leading x — >■ result is obtained by including only the 1-loop vector term, the flrst 
1-loop term in eq. (|3.1|) . The potential is then, equivalently. 



^1-loop/ 63 



-X(j) 
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ISrc'^xy 



(3.3) 
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Two degenerate states are obtained when the last term vanishes. From this one finds 
that 

1 1 

yc{x) = j^-^, y+{x) = —^, y-{x)=0, (3.4) 



^symm 



0, ^broken = ' , (3.5) 

3llX 

1 2^/2 

= 2 2 ' ^3 = 3 5/2 - (3-6) 



Taking into account the other terms in eq. (|3.1|) and the 2-loop effects, one can 
improve on the accuracy. Perturbative results for the critical curve, the latent heat 
and the interface tension in different approximations are shown in Figs. ^ and The 



curves are scaled with the leading a; ^ dependences (eqs. ( |3.4| ), ( |3.6| )). The variation 
as a function of x in the 1-loop results is due to the scalar loop terms in eq. (|3.1| ). The 
2-loop results have been computed by either optimising |T2|, ^ the scale /i or by fixing 
it to be /i = e|. The scale dependence gives an idea of the magnitude of higher order 
corrections. One observes that perturbation theory is becoming unreliable at least for 
x>0.1...0.2. 

Concerning the mass spectrum, using the approximations mw = 630;,, rrijj = V"{(l)b), 
one has in the broken phase. 



mw/ei = l + Jl-167rV /(47rx), (3.7) 



ml 



le\ = (l - ISttV + V^l - Ibix'^yx^ /(Svr^x). (3.8) 



More accurate expressions for these correlator masses as poles of the 1-loop propagators 
are given in p7| . 

The mass spectrum in the symmetric phase is a more subtle issue. The photon 
is argued to be massless to all orders in perturbation theory ||^, Nevertheless, 
non-perturbative effects can produce a mass in the discretized theory, see Sec. 5. The 
scalar excitation, on the other hand, is represented by a bound state in a logarithmically 
confining potential and thus the computation of its mass is non-trivial. A simple way 
to see this is to consider the 3d theory in Minkowskian terms, so that the excitation 
masses are masses of states living in two space dimensions (note that the masses are 
not those of bound states in 3d). The scalar excitation then corresponds to a <^^<^^ 
bound state in a 2d Coulombic potential 

^coui(r) = ^K^{[ir) ~ ^ log — ; (3.9) 
zvr In fir 

/i ~ e| is an infrared (IR) cutoff. A simple minimisation argument then shows that 
the masses are of the form 



M = 2m3 - log ^ + const x el + 0(^), (3.10) 
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Figure 2: The quantity x x yc{x) computed perturbatively for the U(l)+Higgs theory 
to different accuracies. The horizontal hne is the leading a; ^ result in eq. (|3.4| ). The 
lattice Monte Carlo point for x = 0.0463 is also shown. For x = 2 perturbation theory 
is meaningless; then the lattice number is x x ydx = 2) ^ —0.10. 
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Figure 3: The same as Fig. |^ but for the jump of the order parameter (xx^) and for 
the interface tension (xx^/^). 



where the logarithm comes from the Coulombic term and the constant depends on the 
IR cut-off fi and on the quantum numbers of the state. Since the constant is sensitive 
to the IR cut-off, it is to be determined by numerical means. 
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4 U(l)+Higgs theory on the lattice 



Our aim is to study the continuum theory in eq. ( P^ ) and it is thus crucial to keep 
the continuum variables y, x fixed when discretizing the theory. This point has to 
do with ultraviolet renormalisation. The linear and logarithmic mass divergences are 
renormalised in the continuum theory in the MS scheme. Discretization in itself is 
another scheme. The two schemes now have to be related so that correlators measured 
in the lattice scheme go to the correct continuum ones at given y,x m the limit that 
the lattice spacing a goes to zero. This computation has been carried out in [54-56]; 
for order 0{a) corrections, see |p7[ . 

To discretize, we use the compact formulation of the gauge field and introduce the 
link field f/i(x) = exp[zae374j(x)] = exp[iaj(x)]. Rescaling the continuum scalar field 
to a dimensionless lattice field by 



(4.1) 



the lattice action becomes 



S 



1 - cos4 (x)j - /?H E Re 0*(x)f/i(x)0(x + 



X,i<j 



X,i 



+ 0*(x)0(x) +^nY. [</>*(x)0(x) - 1]^ , 



(4.2) 



where [3g = l/e|a, 



aj(x) + aj(x + i) 



Uj (^x } 



(4.3) 



and a specific, customary choice has been made for (3h (other choices are possible as 
well; see for the action before any choice has been made). With these conventions, 
the lattice couplings Pg,Ph,Pr are determined from ^ 





1 

eia 


Pr = 


xf3j, 






2^G 





167r2 



3.1759115(1 + 2x)/?G 
;-4 + 8x- 8x2)(log6/?G + 0.09) + 25.5 + 4.6x 



(4.4) 
(4.5) 



Thus for a given continuum theory depending on one scale 63 and the two dimensionless 
parameters y, x, the use of a lattice introduces a regulator scale a, and eqs. (^^ ) -( ^^) 
specify, up to terms of order CgO the corresponding lattice action. 

In a non-compact formulation of the gauge field one would replace 1 — cos Fij 
in eq. ([4.2|). Then the coefficient 25.5 in eq. (|4.5[) should be replaced with -1.1 |[56||. 



8 



An important quantity we will study is the scalar condensate (0*0). In perturbation 
theory this has a linear and a logarithmic divergence and renormalisation induces a 
known logarithmic scale dependence The condensate itself is not a 

physical quantity, only its changes are. The relation between the continuum (at the 
scale /X = 63) and lattice condensates is 



cont 



/latt 



3.1759115/5g 
An 



log(6/3G) + 0.668 



(4.6) 



where "latt" refers to the normalisation of the field in the lattice action (^4.2| ). The 
first term here is the scaling in eq. (|]1), the second is the linear, and the third the 
logarithmic divergence. The discontinuity of (0*</)(/i)) in a first order transition is free 
from these divergences, as they are just constants. 

For mass measurements we use the discretized forms of the operators discussed in 
Sec. 2. The 0"^"*" and charge conjugation even "scalar" operators (f)*(f) and Re0*Di0 
and the 1 vector operators Im0*Dj0 and eijkFjk are on the lattice represented by 



5(x) = 0*(x)0(x), 

S,(x) = Re0*(x)[/,(x)0(x + ^), 

Vi{^) = Im0*(x)f/i(x)0(x + z), 

Vi(x) = -eij7cSinFjfc(x). 



(4.7) 
(4.8) 
(4.9) 

(4.10) 



In the 4d case these operators have been used for mass measurements in 58, 5£]. Note 



that other (higher dimensional) operators with the same quantum numbers could be 
considered as well, and a systematic way of finding out which combinations couple to 
the lowest physical mass states, is with the blocking and consequent diagonalization 
techniques discussed in Sec. 6.5. In practice, the masses are measured from plane-plane 
correlators instead of the local operators in eqs. (|4.7| )- (|4.1CI|) as discussed in and 



the photon mass measurement requires an external momentum pO|. 



5 Photon in 3d theories 

Various 3d continuum theories with a photon in the physical spectrum are listed in 
Table 1. Here we shall add to this list the U(l)+Higgs theory in the compact lattice 
formulation, eq. ( [4.2|) and, for comparison, also the discretized pure U(l) gauge theory 
in the compact formulation. The situation with photons in these theories is as follows: 
• In 3d SU(2)xU(l)+fundamental Higgs theory there is a massless field, "photon", 
in both the symmetric and broken phases. It is appears both perturbatively and 
non-perturbatively as a physical state [^. In the symmetric phase the photon is 
represented by the hypercharge field, whereas in the broken phase the photon field is a 
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Theory 


d.o.f.'s in 


d.o.f.'s in 


7 in 


7 in 




symm. ph. 


broken ph. 


symm. ph.? 


broken ph.? 


U(l)+Higgs 


A,<f)*,<p 

= 2+2 


= 3+1 


yes 


no 


SU(2)+adj. Higgs 




Wt,A,H 


no 


yes— s-no 




= 3-2+3 


= 2-3+2+1 






SU(2)xU(l)+fund. 




Wr,Zl\A,H 


yes 


yes 


Higgs H 


= 3-2+2+4 


= 3-3+2+1 







Table 1: 3d continuum theories with a massless photon. The 2nd and 3rd columns list 
and count the free field degrees of freedom in the action (with and without a shift to the 
broken minimum) . The last columns express whether a massless photon really appears 
in the physical spectrum. For SU(2)+adjoint Higgs the perturbatively massless photon 
becomes massive through non-perturbative effects [pO[] . 



linear combination of the hypercharge and SU(2) gauge fields. Because the photon is 
massless in both phases, the two phases cannot be distinguished by its value and the 
phase diagram can be smoothly connected. 

• In 3d SU(2)+adjoint Higgs theory there is, perturbatively, a massless photon in 
the broken phase but not in the symmetric phase. This is since after symmetry break- 
ing by (j)"' = (0,0,f) there remains a compact U(l) invariance, related to rotations 
between the 1,2 components. On the tree level the photon could thus be used to 
distinguish the two phases and the phase diagram seems to be split into two disjoint 
regions. However, what happens non-perturbatively [48-51] is that the photon be- 
comes massive due to monopole configurations. The physics of how a gas of monopoles 
induces a mass is familiar from that of Debye screening in a plasma: if there is a gas 
of electrons (e, T, me, ne) in a neutralising background, there is a plasma frequency 
ujp = e^fie/ (meeo) and a screening length proportional to the inverse of Up, so that the 
static photon correlation length has become finite due to charge screening. Since the 
monopole charge is ~ fi/gs, where is the non-abelian gauge coupling, one similarly 
obtains ^ 

~ —^monopole- (5.1) 
93 



One can further carry out a semiclassical computation leading to an expression of 
the type (the exponent is here for x = 0) 



'^monopole exp 



9i 



(5.2) 



Here m^y is the perturbative mass of the the states (Table 1), and one should keep 
in mind that they are actually not physical states (since they are charged). Anyway, 
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the fundamental quantity is the photon mass in eq. (|5.1| ) which can be non-zero and 
has to be measured numerically. This has been done in pT| , |^ and the phase diagram 
has been shown to be smoothly connected. 

• In 3d U(l)+Higgs theory or in pure U(l) gauge theory on the lattice in the com- 
pact formulation there also is a massive photon. A semiclassical computation for this 
Polyakov mass gives in the weak coupling limit /5g ^ 1 5B| 



P I 2 



7r(2/3G)'/'exp 



3.1767r 



G 



(5.3) 



In the strong coupling limit of small (3g-, on the other hand 



P / 2 



(5.4) 



Hence even in the compact U(l)+Higgs case, the phase diagram is connected for all 
— 1/ (3gi and the phase transition only appears in the limit (3g ^ oo. Below 



finite e^a 



we measure explicitly with a series of finite /Sg's and verify this behaviour. Thus, the 
massless state which appears in perturbation theory in the symmetric phase remains 
there in the full non-perturbative case in the continuum limit, as well, dividing the 
phase diagram into two disjoint regions. 



6 Simulations 

In this Section we discuss the numerical results of our lattice simulations. In Sec. 6.1 
we discuss the parameter values used, in Sec. 6.2 the location of the phase transition 
point, in Sec. 6.3 the latent heat of the transition, in Sec. 6.4 our estimates for the 
interface tension, in Sec. 6.5 mass measurements and in Sec. 6.6 critical exponents. 



6.1 Parameter values 

We want to focus our attention on two separate regions: that of small x (eq. ( p.4| )), 
where perturbation theory is expected to work and the transition is of first order, 
and that of large x [Q, where perturbation theory breaks down. In terms of the 
original physical theories these would correspond to type I superconductors or small 
Higgs masses and to type II superconductors or large Higgs masses, respectively. Some 
guidance concerning the limiting value is obtained from the fact that for the other 3d 



theories studied P3| , |21| the first order transition terminates at x = Xc ~ 0.1. ..0.3 and 
that X = 0.5 is the tree- level limiting value between type I and II superconductors. 

Since the study of this problem is rather demanding in computer resources, we have 
chosen for the simulations two values of x, x = 0.0463 and x = 2, corresponding to 
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strongly type I and type II superconductors. When mapped to hot scalar electrody- 
namics, these correspond to rriH = 30, 160 GeV. 

One of the most essential points of the present lattice simulations is that the aim is 
to obtain results for the continuum theory { ^.2\j at fixed y, x. The extrapolation to the 
continuum limit thus has to be carried out carefully: one first takes a fixed a (fixed 
Pg, eq. (|4.5|)) and extrapolates to infinite volume V oo, and then one extrapolates 
a ^ 0. To estimate the required lattice sizes of V = (Na)^ and of the lattice spacing 
a, one notes that the lattices must be big enough to contain all relevant correlation 
lengths, and fine enough so that all the relevant correlation lengths are longer than 
several lattice spacings. Consider a system with a typical correlation length ^. Then, 
accordingly, one has to satisfy (on a periodic lattice) a ^ ^ ^ Na/2 or, in physical 
units, 

ela = 1/Pg « el^ « N/{2Pg). (6.1) 
Thus we have a lower limit for A^: 

N » ej^ ■ 2/3g, (6.2) 

and a lower limit for /Sg'- 

(3g » l/eli. (6.3) 

Eq. (|6.2| ) expresses quantitatively why it is difficult to approach the continuum limit 
/3g ^ oo: the lattice size must be increased simultaneously. We shall use /3g > 4 
and find that typically e\S, varies between 0.3 and 2. Our lattice sizes are thus chosen 
as = 32, . . . , 64. The lattices used are shown in Table ^ 



X 




volumes 


0.0463 


4 


8^ X 40 
16^ X 80^ 


12^ X 60 








6 


8^ X 40 


12^ X 60 










16^ X 80„ 


242 X 90„ 


322 


X 90m 






64^ X 128 










8 


12^ X 60 


16^ X 80 










24^ X 90^ 


32^ X 90„ 


40^ 


X 120m 




12 


16^ X 80 


24^ X 90 


322 


X 90 


2 


4 


32^ 


48^ 


602 


X 120 



Table 2: The lattice sizes used for the simulations at the transition temperature for 
each (x, /3G')-pair. In all cases, several /3iy-values were used around the transition point. 
Multicanonical simulations are marked with the subscript (m)- 

For each lattice listed in Table H we have performed simulations with several values 
of (5h- Typically at x = 0.0463 each lattice has 3-5 different values of jSu and at x = 2 
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somewhat more. For x = 0.0463 where we expect to find a first order transition, the 
different values of f3H are joined together with the Ferrenberg-Swendsen multihistogram 
method |6l|]. This was not done for the x = 2 data. The total number of combined 
heat-bath/overrelaxation sweeps is between 50000 and 150000 for each value. The 
measurements were not performed after every sweep in all cases, so the statistical 
sample sizes range between 20000 and 100000. 

In the first order regime one expects to find supercritical slowing down when one 
goes to large enough lattices. We find that for lattice sizes larger than V/ ~ 300 one 
has to use the multicanonical algorithm to ensure that a correct statistical sample is 
generated. 

The simulations were carried out with a Cray C94 at the Finnish Center for Scientific 
Computing. The total amount of computer power used was about 2500 h of CPU time 
or about 4 x 10^^ fioating point operations ~ 130 Mfiops year. 

6.2 The critical point 
6.2.1 Type I, x = 0.0463 

In principle, if the lattice spacing is small enough, the critical point can be obtained 
directly from Monte Carlo data by observing where the mass of the lowest vector exci- 
tation vanishes. However, in practice the mass measurements are rather demanding 
in the vicinity of the critical point, and if possible we prefer to work with local op- 
erators. There are no known gauge invariant local order parameters, so we use order 
parameter -like quantities, which display a discontinuity at the transition point. The 
actual operators used are B? = (j)*(f) and the hopping term 

L = i^Re 0*(a;)?7,(x)0(x + z), (6.4) 
-J i=i 

averaged over the volume. 

For each individual lattice and (/3g, x) pair we locate the pseudocritical point j3H,c 
using several different methods: 

1. "Equal weight" of the histogram p of R^. 

2. "Equal weight" of the histogram p of L. 

3. Maximum of the susceptibility x of i?^ . 

4. Maximum of the susceptibility x of 

5. Minimum of the 4th order Binder cumulant B oi L. 
We give the susceptibility in dimensionless units as 

= elv {m - m)f) = {{w - mf) , (6.5) 

where the overbar denotes a volume average, = J2x R'^i'^)- The Binder cumu- 
lant B has the standard definition. 
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Figure 4: Equal weight histograms of R'^ and L for x = 0.0463, /3g = 4,6,8,12. As 
the volume is increased the peaks remain at fixed positions as expected for a 1st order 
transition. For Pq — 12 a. proper y — > oo limit would necessitate larger lattices. 
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Figure 5: The susceptibility of divided by the volume V (left) and the Binder 
cumulant of L (right) at x = 0.0463, /3g = 6. 



The equal weight histograms for jSc = 4, 6, 8, 12 are shown in Fig. ^. The histograms 
are obtained by joining all data together and then reweighting it to obtain the point 
where the areas of the two peaks are equal. This involves choosing a value to separate 
the two peaks. The choice is arbitrary, and introduces some systematic error to the 
determination of the critical point. We have used the same value for all lattices for a 
given X, /^c-pair. This value was determined by the minimum of the distribution in the 
largest volume. The error we quote is purely statistical: it was obtained by a jackknife 
analysis. 

The locations of the maxima of the two susceptibilities and the minimum of the 
cumulant were also obtained by joining all data and then reweighting each jackknife 
block independently. As an example, the susceptibility of and the Binder cumulant 
of L for /^G = 6 are plotted in Fig. |]. The maxima of Xr'^ are plotted in Fig. |]. 

In Fig. H we have plotted different pseudocritical points for x = 0.0463, /?g = 4 . . . 12. 
The values obtained from different methods differ at a finite volume, but as one in- 
creases the volume these pseudocritical values come closer to each other and agree with 
each other at the infinite volume limit (within statistical errors). The results obtained 
from different methods are not statistically independent, but serve as a consistency 
check for the infinite volume extrapolation. The different pseudocritical points at the 
infinite volume limit are collected in Table |^. 

The extrapolation of i/c to continuum is shown in Fig. |^. A linear extrapolation has 
a confidence level of 19% and gives the value i/c = 0.1038(14). The quadratic fit gives 
He = 0.120(9). As the linear fit is reasonably good, we choose to quote that value as 
our continuum result. The 2-loop perturbative value t/c = 0.0947 is more than 5a from 
this. 

It is possible to compare results at /5g = 4 and 8 to those obtained with non- 
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x=0.0463, 6^=4 



0.34086 



0.34084 



0.34082 



0.34080 



0.33850 



x=0.0463, P3=6 



• Equal weight of L 
□ Equal weight of 

♦ Maximum of Xi_ 
Maximum of 

M Minimum of B(L) 
Minimum of B(R^) 



4c> 



1/V*e; 
x=0.0463, f 



0.337090 



0.337080 



• Equal weight of L 
□ Equal weight of 

♦ Maximum of Xi 
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M Minimum of B(L) 
, Minimum of B(R^) 
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Figure 6: Different pseudocritical points for x = 0.0463, /3g = 4,6,8, 12 as a function 
of tlie volume. 







piL) 




max XL 


min Bl 


4 


.3408199(18) 


.3408208(19) 


.3408208(25) 


.3408192(28) 


.3408189(34) 


6 


.3383148(07) 


.3383151(06) 


.3383073(37) 


.3383111(24) 


.3383075(23) 


8 


.3370630(10) 


.3370630(10) 


.3370629(16) 


.3370639(23) 


.3370638(23) 


12 


.3358169(40) 


.3358165(40) 


.3358314(107) 


.3358235(92) 


.3358069(86) 



Table 3: The critical coupling (3h,c extrapolated to infinite volume for all lattice spac- 
ings. The first two columns refer to the equal weight criterion. 



compact simulations in ||3^. We find that our critical couplings Ph^c do not agree with 
the critical couplings obtained with noncompact simulations. From Fig. 5 in one 
can read that the critical value for /3g = 4 is roughly 0.34040, with a very small error 
(smaller than 10"^). At /3g = 4 we obtain /3h,c = 0.3408208(19), so the statistical 
error cannot explain this difference. This is not unexpected, as there is no need for 
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0.130. 



0.110 - 



0.100 




Figure 7: The infinite volume extrapolations of the pseudocritical points as a function 
of the lattice spacing ae| = 1//3g- Quadratic and linear fits for the limit a ^ are 
also shown. 



the non-compact and compact results to agree at a finite lattice spacing. However, in 
the continuum limit one should obtain the same results. The authors of ||3^ choose to 
quote the critical temperature instead of the critical coupling Uc. If one extrapolates 
the results from to continuum with a linear fit, one obtains = 130.86(58). Using 
their eq. (13) one can try to compare our continuum value with this non-compact result. 
Our linear fit result corresponds to Tc = 131.28(4), so that the two results agree well 
in the continuum limit. 



6.2.2 Type II: x = 2 



The location of the transition point in the type II regime has already been discussed 
in [38|. Measurements of averages of local operators in the type II regime show no 
evidence of discontinuities. Typical distributions are shown in Fig. To quantify 
the effect, the maxima of the B? susceptibility are shown in Fig. ^ as a function of 
the volume. One can see that the maximum of Xr^ remains constant. This finite size 
scaling behaviour is expected if there is no transition or if the transition is of second 
order with a critical exponent a < 0. For understanding the type II regime and for 
quoting a definite value for ?/c, mass measurements, especially that of m-^, are required. 



6.3 Latent heat 

The latent heat — or the energy released in the transition 
to M-i = A (0*0) (see, e.g.. Sec. 11 in |2^ 
perturbatively from lattice simulations. In Fig. |10| we have plotted the infinite volume 



is directly proportional 
This quantity is easily computed non- 
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Figure 8: Histograms for i?^ at x = 2, around the critical point y 
two-peak structure is seen. 
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Figure 9: The maximum of the susceptibility of for x = 0.0463 and x 
function of the volume. 
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extrapolations of A (i?^). These can be converted to continuum values by eq. ([4.1|), 

A (0*0) 1 



(6.6) 



The infinite volume limit is taken by extrapolating linearly with the inverse area of the 
system, as extrapolation with the inverse volume would fail to accommodate the small 
differences between the different ratios of N^/Nz- 



The continuum limit of Ai^ is displayed in Fig. The continuum extrapolation is 
done with a linear fit, which seems to fit the data well - the confidence level of this fit 
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is 28%. The final continuum value we obtain is A^s/cg = 1.64(12), with a statistical 
error only. 

Again the results should be compared with those obtained from perturbation theory 



and from non-compact lattice simulations. In ||3^ it was found that the value of A^s 
increases with decreasing lattice spacing; however we find a rather strong decrease. 
Therefore, even though all our data except jSc = 12 give a value higher than the 
perturbative one, our continuum value is considerably smaller than the perturbative 
value A£3 = 2.25. Taking into account the fact that our data at Pg = 12 is of not 
as high quality as at smaller lattice spacings by using only the three largest lattice 
spacings for the continuum extrapolation, one still obtains a value smaller than the 
perturbative one, Ai^, = 1.71(13). 

6.4 Interface tension 

The interface tension is the most difficult of the characteristics of the first order phase 
transition to measure. We estimate it with the histogram method at the critical 
point the system can reside in a mixed state consisting of domains of the pure phases. 
Because some extra free energy is associated with the interface separating the two 
phases, the area of the interface tends to a minimum and in a system with a geometry 
A^^ X Nz (with A^^ > A"^), the minimum area is just 2 x N^a^ (due to periodic boundary 
conditions at least two interfaces are needed). Since the cost of an interface is ~ 
exp[—aA/T], the dimensionless interface tension can be calculated from 



1 P 

, max ( „ „N 
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Figure 11: The infinite volume limit (left) and the extrapolated results at different 
lattice spacings (right) for the dimensionless interface tension a^. The systematic 
errors are expected to be much larger than the statistical ones shown in the figures. 



where Pmin and Pmax are the minimum and the maximum of the probability distribu- 
tion, and 



eiA 



eiN-'a 



2^2 



'G . 



For finite size corrections in the interface tension, see, e.g., p 



(6.8) 



The formula ( |6.7|) requires that the two interfaces are far enough from each other 
so that their mutual interaction can be neglected. This would be signalled by a flat 
minimum in order parameter histograms, which should appear at large enough values 
of Nz- As seen from Fig. ^, a clearly flat minimum is not yet observed: the interfaces 
are so thick that their interaction is not negligible even for lattices of length 90. Thus 
non-negligible systematic errors are expected in our results. 

The infinite volume extrapolations and the values obtained at different lattice spac- 
ings are shown in Figure |Tl|. No extrapolation to continuum limit is performed even 
though the data seems to be consistent with a linear fit. This is since a linear fit 
would give a negative value for as, which is unacceptable. We ascribe this unphysical 
behaviour to larger systematic errors at larger /^c's, where the physical volume is not 
large enough to contain regions of clearly separated phases. The results we obtain are 
as = 0.28(2) at (3g = 4, cxg = 0.15(1) at I3g = 6, and = 0.09(2) at (3g = 8; (3g = 12 
does not give a good enough signal for this purpose. 

The perturbative value for cxs is 0.225, and as the results at /9g = 6 and 8 are both 
below this and the value seems to be decreasing with decreasing lattice spacing, one 
expects that the real value is significantly lower than the perturbative one. In Fig. |^ 
we have shown the conservative estimate cxs = 0.14(14) which contains all the values 
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Figure 12: Effective masses at different bloclcing levels at x = 2, /^c = 4 as a function 
of z (in lattice units). 



we have measured. 



6.5 Mass measurements 

Since bulk quantities are rather insensitive to the phase transition at large values of x, 
we have paid special attention to mass (= inverse correlation length) measurements. 
Most of the results of our measurements have been reported in [Q, and here we 
describe the techniques used in some more detail. There are several factors which 
could result in too high a value for a mass. First, one has to be sure that there is 
no contamination from higher excited states. This can be ensured by systematically 
searching for optimized operators coupling to the desired excitations, using blocking 
and diagonalization techniques. One also has to start the fits from sufficiently large 
values of z and to monitor the behaviour of the effective mass. Second, one should use 
large enough lattices — in particular, one cannot expect to measure reliably masses 
smaller than ~ 1/iV. 



6.5.1 Blocking 



To avoid having to use extremely large starting values of z it is necessary to have a 
very good projection on the ground state. The projection can be improved upon by 
using blocking in the transverse direction |2^, |6^ both for the link variables Ui{x), 



■n+l. 



X 



u^xwnx + + V uuxwnx + j)unx + j + twux + 2ty 



(6.9) 
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and for the scalar fields 0(x), 

1 r ^ .1 

(6.10) 



i=±i 

The operators measured are then constructed from the blocked fields according to 
eqs. ( |4.7| )-( P?T0| ). The iteration level n can be tuned for optimal overlap. In Fig. |12 



we show the effect of blocking. Typically we find that the best results are obtained at 
blocking level 3 for all operators. The overlap of our blocked operators with the ground 
state ranged from being consistent with 100% far away from the critical point to 75% 
at the critical point. 

6.5.2 Variational analysis 

The blocking described above is in practice sufficient for the scalar excitations, where 
we are mostly interested in the ground state. However, we would like to know not 
only the ground state of the vector excitation but also have some information on the 
excited state. We are especially interested in knowing to what extent the lowest vector 
excitation couples to the operator which consists purely of gauge fields. To study 
this we have performed a variational analysis, which enables us to measure also the 
excited state. 

Variational analysis is based on the fact that we have several operators which couple 
to the same quantum numbers. Even if each individual operator is not perfect in 
the sense that it also contains contributions from higher excited states, it should be 
possible to reduce this effect by suitably choosing linear combinations of individual 
operators. The basis that we use consists of different operators measured at different 
blocking levels. To be specific we use blocking levels 2 and 3; for scalar excitations we 
use the operators 0*0 and Re0*f/j0 and for vector excitations Im0*f/j0 and sin Fjk, 
see eqs. (|47| )- (|4.10| ). Thus in both cases we have four different operators, which we 



call Oj. We then measure the whole 4x4 cross correlation matrix 

a,{z) = {o:{z)o,m, (e.n) 

and then construct the improved operators Oi in the standard way [p5|], 

0, = J2vlO,. (6.12) 

i 

The operators are further normalised to unity at zero distance. The coefficients of the 
operator Oi for the ground state are found by maximising 

(Oi*(a)Oi(0)). (6.13) 

The operators for excited states can then be obtained by maximising the analogy of 
eq. ( |6.13| ) in a subspace that is orthogonal to the ground state operator Oi. 

We find that the variational analysis has usually little effect for the ground state — 
however, for excited states it is extremely helpful. 
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6.5.3 Fits 



We extract the different masses by fitting exponentials to tlie corresponding correlation 
functions. In practice we have used one and two exponential fits and find that the 
difference between these two is within statistical errors. However, the one exponential 
fit is usually more stable, and therefore all numbers we quote are from one exponential 
fits. 

One should note that the expectation values of the correlation function C{z) are 
strongly correlated for close values of z. Typically the statistical correlation between 
C{z) and C{z + 1) is greater than 95%. Therefore it is necessary to take this correlation 
into account when performing the fits. We find that the effect of taking this into account 
is to constrain the fits more tightly. The is increased and the variation between 
different jackknife blocks is decreased. If we do not use correlation information in our 
fits the fit parameters can fluctuate more freely causing not only unrealistically large 
errors but also changing the final value of the fit. 

The major systematic error comes from the choice of the fitting range. We have 
carefully monitored the behaviour of the local mass, and start the fits only when a 
clear plateau is visible. Still changing the starting point by one or two lattice spacings 
changes the result; however as the change typically is roughly equal to the statistical 
error, we do not quote a systematic error on our fits. When interpreting the results 
one should bear in mind that it is possible that the systematic errors are at least as 
large as the statistical ones. 



6.5.4 Polyakov mass 

Let us apply the methods described first to the photon mass in the compact formulation, 
see eqs. (^.3|), (|5l^). In order to use as an order parameter, one has to be sure that 
this Polyakov mass it acquires at a finite lattice spacing, is negligible for all practical 
purposes. The semiclassical calculation ( p. 3D gives m^/e\ = 0.0033 at Pg = 4, which 
should be much smaller than anything we expect to be able to measure. However, as 
most of our conclusions depend strongly on m^, we have checked the vahdity of the 
semiclassical calculation by direct simulations in the pure U(l) gauge theory. We used 
a 24^ lattice, varied /5g = 1 . . . 4 and used the operator (^4.10|) to measure m^. 

The Monte Carlo results for pure gauge theory are shown in Fig. |T3|, together with 
the semiclassical result ( ^.3|) and three data points in the U(l)+Higgs theory. The 
U(l)+Higgs data points were taken in the symmetric phase: we chose x = 2,y = 
and y = 0.5. One sees that adding the Higgs field does not affect the Polyakov 
mass, and that the Monte Carlo data agree with the semiclassical result for /?g ^ 3 
(in agreement with semiclassical computations was found even at Pg ~ 2, but 
this relies on an "improved" Pg- Using that, the agreement in Fig. ^ is improved, 
as well). At /9g = 4 the mass is clearly too small to be seen, and therefore can be 
neglected for all practical purposes. For lower values of Pg we find that the formula 
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Figure 13: The photon mass as a function of Pc in compact 3d U(l) gauge theory. The 
5 X semiclassical curve extrapolates between weak and strong coupling regimes, and 
reproduces the non-perturbative values reasonably well for intermediate Pc- 
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Figure 14: Masses near the critical point at x = 0.0463, /^g = 6 (left) and x = 2, (3g = 4 
(right) (from p8| ). 



5 X semiclassical result reproduces the Monte Carlo data reasonably well, effectively 
interpolating between the strong and weak coupling results in eqs. ( |5.3| ), ( |5.4] ). 

6.5.5 Masses 

The masses in the U(l)+Higgs theory obtained with the methods described above are 
plotted in Fig. |T^. We see that the mass of the photon vanishes within 2a for both 
X = 0.0463 and x = 2 in the symmetric phase. When measured in ex's, some of the 
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results at X = 0.0463 are somewhat further from zero. We beheve that this due to the 
fact that no variational analysis was done for these points. Thus serves as an order 
parameter for practical purposes. 

The first order nature of the transition at x = 0.0463 is clearly visible. All the masses 
show a clear discontinuity at ydx). When going from the broken to the symmetric 
phase, rriH jumps upwards (the scalar correlation length decreases from about 12a to 
about 7a). The two vector excitations and mw, which were degenerate for y < yc, 
separate so that jumps to zero and mw becomes the mass of an excited two-scalar 
bound state. 

At X = 2 the pattern is less clear. The mass of the scalar excitation tuh dips deeply, 
but does not go to zero. The value of rriH at the minimum was also shown to be 



independent of the volume within statistical errors In the broken phase the two 
vector operators again see one single excitation, but in the symmetric phase has 
dropped to zero and the other excitation has a large mass. However, on the basis of this 
data one cannot conclusively state what is happening. One can envisage the following 
possibihties: 

• There is one line ydx) on which goes continuously to zero, mn has a minimum 
and m^,mw separate. 

• As the previous one but jumps discontinuously to zero. This does not imply 
that the transition is of first order, as there is no jump in local quantities. 

• The line ydx) has split into several branches. 

• All the masses go to zero at ydx). This is the standard 2nd order scenario, which 
however, is not supported by our data. 

The data shown is at fixed /^g = 4 and one may ask whether the extrapolation to 
continuum would change the conclusions. The minimum of niH in Fig. corresponds 
to a correlation length of 16a and it is hard to see how making a still smaller could 
change the results. We have made some runs at /5g = 3, 6 and did indeed not find any 
dependence outside of statistical errors. 



6.6 Critical exponents 

If goes to zero continuously, one can measure its critical exponents. If the transition 
were of second order in the usual sense, there would be critical exponents for both scalar 
and vector correlation lengths: 

mn ~ \y-ycr, (6.14) 
~ \y — yc\'^. (6.15) 
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y range 


A 




V 


CL 


-0.25... 


1.07(04) 


-0.040(05) 


0.5* 


0.14 


-0.25... 


0.91(24) 


-0.055(23) 


0.39(17) 


0.10 


-0.5... 0.5 


0.995(19) 


-0.033(04) 


0.5* 


0.08 


-0.5... 0.5 


0.931(62) 


-0.048(13) 


0.429(62) 


0.07 


-1.25... 0.5 


0.959(11) 


-0.026(04) 


0.5* 


0.006 


-1.25... 0.5 


0.947(12) 


-0.046(06) 


0.440(19) 


0.06 



Table 4: Fits to determine the critical exponent v. The superscript (*) after a value 
means that it was kept fixed. 



However, as we do not see any critical behaviour for the scalar correlation length, we 
cannot measure any critical exponent corresponding to it. The approximate scaling of 



mil away from the critical point has been discussed in ||38|| . 

To measure the critical exponent v corresponding to a diverging photon correlation 
length, we have fitted functions of the form vrir^je^ = A{yc — yY to the data shown 



in Fig. 14, both keeping the critical exponent y fixed to the mean field value 0.5 or 



allowing it to vary freely. We have also varied the range of y which we include in 
the fits as ymin = —1.25. . . — 0.25. The fits with their confidence levels are shown 
in Table |[ We find that the critical point is rather stable, and that the value of the 
critical exponent is consistent with the mean field value. All the fits have rather low 
confidence levels, but as noted previously the errors we quote for masses are statistical 
only. Inclusion of the systematical errors (which we estimate to be of the same order 
of magnitude as the statistical errors) should bring the confidence levels to perfectly 
acceptable values. 

In principle one could also measure the anomalous dimension rj^ oi the photon 
correlation length, (V3(x)V3(0)) ~ l/|x|^+''^, but this would require a significantly 



extended analysis and more accurate data. 



7 Conclusions 

We have in this paper studied the 3d U(l)+Higgs model as an effective theory of 
various finite temperature theories. There are two separate aspects of the problem: 
computing the parameters of the effective theory in terms of the physical parameters 
of the full theory and studying the 3d theory as such. The first aspect has been solved 
for a very concrete physical phenomenon, superconductivity (Appendix A), and for a 
more academic system, hot scalar electrodynamics (Appendix B). 

The second aspect, the study of the 3d theory as such, also splits into two parts: the 
region of first order transitions, corresponding to small x = Aa/cg (type I supercon- 
ductors, or small Higgs masses), and the region of possibly second order transitions. 
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large x (type II superconductors, or large > mw Higgs masses). The type II regime 
has been analysed in |^8|], and here we have mainly concentrated on the type I regime. 

The evidence for a first order phase transition at x = 0.0463 is compelling: several 
order parameter like quantities display a clear discontinuity, the finite volume scaling of 
the susceptibities is clearly consistent with a first order transition and the mass of the 
photon can be used as an effective order parameter. However, we have found that even 
when the transition is relatively strong, there are quantitatively some discrepancies 
with perturbation theory which seems to predict too strong a transition, see Fig. ^. 
This is in contrast to the case of the SU(2)+fundamental Higgs model, for example, 
where 2-loop perturbation theory works quite well for comparative Higgs masses |^ . 



The reason is probably just that the transition is weaker in U(l)+Higgs, but one might 
also speculate that physical topological defects (vortices) may play a role. In any case, 
the transition is too weak to be observed to be of first order in practical superconductor 
experiments. 

At X = 2 the situation is even more subtle. Many of the operators that can be 
used to study the phase transition at x = 0.0463 do not display any indication of a 
transition at x = 2. None of the local parameters [R"^ and L) show any discontinuity. 
The maxima of the susceptibilities remain constant. The scalar correlation length 
increases, but still seems to remain finite when extrapolated to infinite volume and 



zero lattice spacing |Q . The only indication of the transition is a vector correlation 
length diverging within statistical errors. 

One thing to be checked in this argument is the use of the compact formulation 
which, in fact, makes the photon massive at a finite lattice spacing, if only by an 
exponentially small amount. Thus in a strict sense the transition vanishes altogether 
for a finite a. We have therefore studied the topological effects related to the compact 
formulation and demonstrated that they behave as expected, vanishing rapidly when 
approaching the continuum limit. Thus we believe that they do not alter the pattern 
described above. 

The precise properties of the transition at x = 2 remain unclear. Our data cannot 
distinguish between a second order phase transition in which the mass of the photon 
vanishes continuously and a more exotic scenario in which the photon mass has a 
discontinuity. To solve this unambiguously one would clearly need much more data. 
Another study requiring even more data would be that of the endpoint Xc of the first 
order regime. 

Because of the non-abelian gauge structure one might have expected computations in 
3d SU(2)+Higgs theories [10,22-28] to be more complicated than those in U(l)+Higgs 
theories. However, just the opposite is the case. U(l)+Higgs is more demanding to 
simulate than SU(2)+Higgs at least in the compact formulation, as some of the corre- 
lation lengths are larger, the transition is weaker, the structure of the phase diagram is 
more complicated, and very large lattices are needed. This is analoguous, say, to what 
happens in the g-state Potts model in which the strength of the transition rapidly in- 
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creases with the number q of field components. Perhaps also the formation of physical 
topological defects, viz. vortices, plays a role in U(l)+Higgs. 
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Appendix A: Full theory I: BCS superconductiv- 
ity and the Ginzburg-Landau model 

So far we have only treated the 3d effective theory as such; in these two Appendices 
we shall quantitatively discuss how two different physical theories map to the same 
effective theory. The first case is superconductivity. 

Quantum phenomena in superconductors are microscopically described by the BCS 
theory In the normal state electrons do not form bound coherent states because 
of the repulsive Coulomb force but at low temperatures the electrons can form Cooper 
pairs due to their interaction with the ionic lattice. The electron-phonon interaction is 
important only at temperatures where the thermal excitations from the Fermi energy 
of the electron gas have an energy smaller than the average phonon energy. When this 
condition is satisfied the Fermi sea becomes unstable against the formation of bound 
pairs of electrons from states above the Fermi surface. 

When the BCS theory is treated within the framework of the Green's functions 
method one can find the connection between the microscopic BCS theory and the 
effective macroscopic Ginzburg-Landau model of superconductivity. Using the electron- 
electron interaction of the BCS theory and solving the equations of motion of the 
Hamiltonian Gor'kov [^] derived an equation describing the free energy density of 
Cooper pairs (we keep here c, %): 
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where A(x) = ('?/;-f(x)'?/'j^(x)), r < 1 (= 
of the material, 



-1 for perfect crystals) is an impurity parameter 
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(A.2) 
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sets the typical length scale of the system, and 



iV(0) 



2^ 



(A.3) 



is the density of states at the Fermi surface. Note that the temperature parameter Tq 
equals only on the mean field level. The connection between the free energy density 
( |A.lj ) and the parameters in eq. ( p.4|) of the U(l)+Higgs theory can now be established 
by scaling the fields so that ( |A.1| ) has the form in eq. ( ^^.21 ). The final relation becomes 
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with Tp = mvp/2 ^ T^. What now is crucial is that widely different scales with large 
or small ratios appear. For low temperature superconductors (Tc ~ 1 K), T^/Tp ~ 10~^ 



and the Fermi velocity is vp ~ (10 



10" 



c. 



These lead to values (7 ~ 10 and 



q ~ 0.01 resulting in x ~ O.Ol/r^. For usual superconductors thus x <^ 1. In contrast, 
for high temperature superconductors typically x ^ 1. This follows from measurements 
of the penetration depth and the coherence length and from the fact that x is simply 
half the square of their ratio, so that x < 1/2 describes type I superconductors and 
X > 1/2 type II superconductors. It is interesting that the physically challenging case 
(high Tc superconductors) just corresponds to the region in which the G-L model has a 
particularly subtle structure. It is so in the relativistic case (Appendix B), as well, that 
the large x regime correponds to the most challenging domain of the physical theory: 
large Higgs masses. 



The big dimensionless ratios also appear in the expression for y in eq. ( A.4 ). This 
implies that although ydx) is of order one for small x, T^/Tq — 1 is extremely small, 
of the order of ~ 10^^. This smallness follows not from the dynamics of the 3d 
theory, but from the relation of the full and effective theories. 

Similar statements apply in the case of the latent heat L. In the effective theory, 
L is essentially the jump Ai^ in (0*0) at the phase transition point. For small x (see 
eq. (|3^); now again h = c = 1), 
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The conversion into physical units gives then 



L = Ai.elT!^ = T^-pi-^ ^ 1.15 x 10-«ir,^-^J/cm^ (A.7) 

where Cp = 2e = a/IStto ~ 0.61 is the couphng constant of the physical theory and Kc 
is the transition temperature in Kelvins. For example, for Pb [K^ = 7.2 K) the latent 
heat is 

L ~ lO^'^J/cm^ = 10~^J/mole, (A.8) 

which is far too small to be detected experimentally at present. This is so in spite 
of the fact that the transition is quite strong from the point of view of the effective 
theory. Note also that non-perturbatively the transition is somewhat weaker than given 
by eq. (|A.6|) , see Fig. |, but this does not change the order of magnitude estimate in 
eq. 



Appendix B: Full theory II: hot relativistic scalar 

+ fermion electrodynamics 

Secondly, we shall work out the 3d effective theory of finite T relativistic scalar electro- 
dynamics. Since we aim at precision, we here first have to discuss how the parameters 
of the Lagrangian are determined to 1-loop in the MS scheme in terms of the mea- 



surable (in principle!) physical masses of the theory |13|. Next we have to relate this 
4d finite T theory and the 3d effective theory. This happens in two stages: first the 
momenta p> vrT (the superheavy modes) are integrated out, resulting in an effective 
bosonic 3d theory with fundamental and adjoint Higgs fields. Finally, the adjoint Higgs 
field (the heavy modes) is integrated out. Some parts of the latter two steps have been 



already discussed in W% 31 1. The parametric accuracies of the effective theories 



derived have been discussed in |T^. 



B.l Relation between physics and the Lagrangian in MS in 4d 

The Lagrangian of 4d scalar + chiral fermion electrodynamics is 

£ = ^F^,F^, + (D^0)*(D^0)-z/V0 + A(0»' 

+^iP^L + ^R0R + 9Y^{(j)*aL + (pan)^, (B.l) 



where 



^H,L = aH,L^, a^„i = (l±75)/2, = (0i + ^02)/V2. (B.2) 
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At the tree-level, the couplings u'^^X^gy are related to the physical masses m//, mw 
and nit of the 0i, A^, and \l/ fields in the broken phase by 
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[B.3) 



The masses mn, mw, are analogous to the Higgs, W and top masses of the Standard 
Model, and they can, in principle, be experimentally determined. The gauge coupling 
is, in principle, determined in terms of some cross section. 

When one takes into account loop corrections, the theory requires renormalization, 
and the relations in eq. ( |B.3| ) change. Let us for definiteness regularize the theory in 
the MS scheme in ci = 4 — 2e dimensions. At 1-loop level, the bare quantites are then 
related to the running renormalized quantities through 



A 



B 



B 



^B 



V 



B 



(f){n) 
v]/(/i) 



-2e 



2e 



1 + 



1 + 



/i 



2e 



(- 

(-2 



,Ns + 2N 



■3e2 _ 9l^ 



6 

2'r 



/i 



-2e 



91b 



A(/i) + 



fx 



'2t 



4 ^ 



-3e2 + 4A + c/^ 



1 + 



(3e^ - 6Xe^ + lOA^ + 2A^^ - 
(2gl - 3e' 



(B.4) 
(B.5) 
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(B.8) 
(B.9) 
(B.IO) 



Here A^s = 1 is the number of scalar fields and Nf = 1 is the number of fermion 
fields; Ns and Nf are used just to show the origin of the different contributions. The 
renormalized parameters of the theory run at 1-loop level as 
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To relate the running parameters to the physical observables, one has to calculate 
suitable physical quantities to 1-loop order, and express the running parameters in 
terms of these. For illustration, let us consider relating the running parameter z/^(/i) 
to physical masses. To do so, we calculate the physical pole mass niH to 1-loop order. 
For the calculation, the field 0i is shifted to the classical broken minimum 0^ = '^^/A. 
Then the radiatively corrected renormalized 1-loop propagator of the Higgs field is of 
the form 

where ml = 2z/^. Solving for the pole niH, one gets 
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1 + 
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(B.16) 



This expression is gauge-independent, since the self-energy is evaluated at the pole. The 
latter term in eq. ( p.l6|) is the 1-loop correction to the tree-level formula in eq. (p.3|) , 
and contains the //-dependence in eq. ( |B.12p . 

To give the explicit form of the 1-loop expression for we use the notation 



h 



mw ' 
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mw 



[B.17) 



Inside 1-loop formulas, one may then write A = e^/i^/2, Qy = 2e^t^. Let us also 
introduce the function 



= -2^4/1^2 _ I arctan 
which has the special value Li{l) - 



1 



/Vs. Then we get 



[B.18) 
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B.2 Integration over the superheavy scale 

At high temperatures, the equilibrium thermodynamics of the theory defined by 



eq. (|B.1|) can be described by a superrenormalizable 3d effective field theory (for a 
review, see [^). The Lagrangian of the effective theory is that in eq. (|2.2| ) with a 
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second scalar Aq. 



£eff = -i^,,F,, + (A0)*(A0)+m20*0 + A3(0*0)' 
+^{d,Aoy + \mlAl + -\aAI + h;<\>^Al 



(B.20) 



where all the fields have the dimension GeV^^^ and the couplings e|, A3, A^, /is have the 
dimension GeV. 

The purpose of this section is to give the parameters 63, /is, mg, A3, m|,, A^ and the 
fields Aj, in terms of the parameters in eq. ( |B.1| ), with the accuracy compatible with 
1-loop renormalization of the vacuum theory, as described above. Apart from terms 
arising from fermions, the results have been given at 1-loop order in (2-loop results 
for the mass parameters have been added in |31[]). Here we add the effects of fermions. 
We use extensively the generic results from |jl3|. The basic notations are: 
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47rT 7.0555r' ' ' ttT 1.7639r 

Let us start with the relations of the wave functions. Using the Z'^^s from eqs. (35- 
37) of H^, the momentum-dependent contribution of the superheavy modes to the 
two-point scalar correlator is 

1 



For the spatial and temporal components of the gauge fields one gets 

Z = Zfjfj^ e NgZ^i^ e N^Z-p-pi 
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where the Z^^s are from eqs. (38-45) of [|1^]. Hence the wave functions in the 3d action 
are related to the renormalized 4d wave functions in the MS scheme by 
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, (B.25) 
(B.26) 



The couplings Cg and /is can be extracted from the superheavy contributions to the 



Aj)- and (0074oAo)-correlators at vanishing external momenta, respectively. The 
contributions from the relevant diagrams are in eqs. (50-63) of [|T^. The result is 

G"" = + Q^^ - ^Xe^Qh - ^e^Giv + ^Xe^Giss ' e^alG^FFF- (B.27) 
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When the fields are redefined according to eqs. ( [B.24|) - (|B.26| ) and the vertex is identified 
with the corresponding vertex in eq. ( |B.20| ), one gets 
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For the scalar sector, the required correlators are most conveniently generated from 
the effective potential. With the Lagrangian masses 
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the 1-loop (unresummed) contribution to the effective potential in Landau gauge is 

V^i(0) = Cs(mi) + Cs{ni2) + Cv(mr) + Cpinij), (B.31) 

where the C's are from eqs. (69-71) of |Il3[. From the term quartic in masses in Vi(0), 
one gets the superheavy contribution to the four-point scalar correlator at vanishing 
momenta. Redefining the field according to eq. ( |B.24| ), the coupling constant A3 is 



then 
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(B.32) 



The coefficient of in Vi(0) gives the 1-loop result for the scalar mass squared. The 
result is the term of order g'^ on the first line of eq. ( |B.36| ). 

For the 2-loop contribution to the mass squared one needs the 2-loop effective 
potential V2(0). In terms of eqs. (81-93) of |]T3|, the result is 
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With the abbreviations 
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the result for the 3d mass parameter is 
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Here we have taken into account higher-order corrections on the last line, by replacing 
the 4d coupling constants by the 3d ones which appear in the exact running of 
inside the superrenormalizable 3d theory. 

The parameters and require the calculation of the superheavy contributions 
on the (AqAq)- and (y4ov4oAoAo)-correlators at vanishing momenta. Using eqs. (96-101) 
of the 1-loop 2-point correlator for the U(l)-field is 



and using eqs. (102-109) of |T3[, the four-point correlator is 
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Since there is no tree-level term corresponding to the correlators in eqs. (p.37| )- (|B.38|) , 
the redefinition of fields in eq. (|B.25|) produces terms of higher order. The final results 
can then be read directly from eqs. (p.37| )- ([B.38|) : 

'''^'(2iV, + iV;), (B.39) 
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Apart from fermionic contributions, the 2-loop result for mf) has been given in |3l 

Using eqs. ([B.11| )- (|B.14D , one sees that the quantities e|, h^, A3, i/^ and gy are 
independent of /i to the order they are presented above. In other words, when the 
running parameters e^(yu), z^^(/u), A(/i) and gy^fx) are expressed in terms of physical 
parameters, the //-dependence cancels in the 3d parameters. Note also that the bare 
mass parameter mg of the 3d theory is independent of /i, so that one may use an 
independent renormalization scale inside the 3d theory. 
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B.3 Integration over the heavy scale 



Consider now the parametric magnitude of the couphngs in eq. ( [B.20|) . If we are 
interested in the study of the phase transition, we know this happens when m^{fi) is 
very small and we see from eq. (|B.36|) that the tree-level term and the 1-loop term 

e^T^. At the same time the mass of Ao is of the 



have to cancel so as to leave m| 



order of eT. The action in eq. ( p.2CI|) can then be further simplified by integrating out 
the Ao-field, leaving an action of the G-L form in eq. ( |2.2| ). 

If we proceed from the T ~ Tc region to T ^ Tc, eq. ( B.36 ) implies that also 
is of the order of eT and thus it should be integrated out together with Aq. The 
resulting theory is very simple: free electrodynamics in 3d. This is the statement that 
the magnetic sector of hot scalar ED is trivial. 

For the case T ~ when the action is the one in eq. (|2.2| ), we denote the new 
parameters by a bar. The relations of the old and the barred parameters have been 
and we include the results here just for completeness. 
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given m 

The integration over the heavy scale proceeds in complete analogy with integration 
over the superheavy scale in Sec. [B.2| . Since there are no momentum-dependent 1-loop 
corrections to the 03- and Aj-fields from the Ao-field, 03 = 03 and Ai = A'f'^. Since 
Aq and Ai do not interact, there are no 1-loop corrections to the (03)^(y4?°')^-vertex, so 
that 63 = 63. The scalar self-couplings get modified, and at 1-loop order 
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To calculate the 2-loop corrections to the mass parameter, we use again the effective 
potential. The 2-loop contribution from the heavy scale to the effective potential is 



(B.42) 



where the D's are from eqs. (119-122) of [|l^. The mass parameter rnKfi) is then 
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where we used eq. (|B.36| ) and included higher order corrections on the last line, by 
replacing the couplings with those of the effective theory. 
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Using the results of this section for and A3 and adding the running parameters 

from Sec. [B.l| , the final parameters y and x of the G-L action in eq. ( |2.2| ) are completely 



fixed in terms of the physical quantities e^, the temperature T, and the physical pole 
masses mw,mH and rrit at zero temperature. We do not write out the expressions 
explicitly. For illustration, let us show the final expression for u"^ from eqs. ( |B.19|) , 
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